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^ : ABSTRACT 

We describe new, deep MIPS photometry and new high signal-to-noise optical 
^ I spectroscopy of the 2.5 Myr-old IC 348 Nebula. To probe the properties of the 

IC 348 disk population, we combine these data with previous optical/infrared 
r-| I photometry and spectroscopy to identify stars with gas accretion, to examine 

^1 their mid-IR colors, and to model their spectral energy distributions. IC 348 

O I contains many sources in different evolutionary states, including protostars and 

^ ' stars surrounded by primordial disks, two kinds of transitional disks, and debris 



OO 
O 



X 



disks. Most disks surrounding eary /intermediate spectral-type stars (> 1.4 M© 
at 2.5 Myr) are debris disks; most disks surrounding solar and subsolar-mass stars 



J> I are primordial disks. At the 1-2 a level, more massive stars also have a smaller 

^ ' frequency of gas accretion and smaller mid-IR luminosities than lower-mass stars. 

These trends are suggestive of a stellar mass- dependent evolution of disks, where 
most disks around high/intermediate- mass stars shed their primordial disks on 
O ' rapid, 2.5 Myr timescales. The frequency of MlPS-detected transitional disks is 

~ 15-35% for stars plausibly more massive than 0.5 M©. The relative frequency 
of transitional disks in IC 348 compared to that for 1 Myr-old Taurus and 5 
Myr-old NGC 2362 is consistent with a transition timescale that is a significant 



' fraction of the total primordial disk lifetime. 

Subject headings: stars: pre-main-sequence-planetary systems: formation 
planetary systems: protoplanetary disks 



1. Introduction 

Most < 1 Myr-old stars are surrounded by optically-thick, accr eting primordial disks, 



which produce strong near-to-mid infrared (IR) emission (L^^/L^ > 0.1: lKenyon and Hartmann 
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19951 ). The gas and dust in these disks comprise the building blocks of planets. As stars 
age, the dust grains grow, settle toward the midplane, and bec ome incorporated into plan- 
etesimals. Circumstellar gas depletes by accretion onto the star (IHartmann et al.lll998l ). All 
of these processes re duce the amount and frequency of disk ernission on tinaescale s of ~ 3-7 
Myr for most stars (IHernandez et al. 1120071 : iHillenbrandl l2008l : ICurrie et al.l 120091 ). 



B y ~ 10 Myr, mos t stars do not have optically-thick primordial disks (ICurrie et al. 



2007a 



Hillenbrandll2008l ) . Disks around these older stars are typically debris disks (ICurrie et al 
2008al ). which have weaker, optically-thin dust emission (L^/L^ < 10"'^) and lack evidence 



for circumstellar gas accretion. Because stellar radiation - radiation pressure and Poynting- 
Robertson drag - or stellar wind drag can remove dust in debris disks on timescales much 
less than the age of the star, their dust re quires a replenishment source, which is supplied by 



activ e icy o r rocky planet formation (e.g. iBackman and Parescelll993l : iKenyon and Bromley 



2004 . 12OO8I ). Because debris disks lack copious amounts of circumstellar gas needed to form 
gas giant planets, identifying when most primordial disks turn into debris disks pinpoints an 
empirical upper limit for the formation timescale for gas giant planets. 

Recent Spitzer Space Telescope studies of 1-5 Myr-old clusters indicate that primordial 
disks disappear faster arou nd early type, high/iii t ermediate-mass s t ars than they do around 



late type, low-ma ss stars (ICarpenter et al.l l2006l : iLada et al.l l2006l : IHernandez et al. 1 12007 



Currie et al.ll2009l ). By 5 Myr, debris disks completely dominate the disk population around 
high/intermediate- mass stars. At 5 Myr, most disk s around lower-mass stars appear to have 



inner holes and/or extremely low warm dust masses (iDahm and Carpenterll2009l : ICurrie et al. 
20091), which indicates that they may be disappearing and actively making the primordial- 
to-debris disk transition. Thus, the epoch at which most high/intermediate- mass stars make 
the primordial-to-debris disk transition must occur at an age prior to 5 Myr and may be 
identified by mid-IR observations of younger clusters. 



With a median age of ~ 2-3 Myr, stars in the IC 348 Nebula (lHerbig||l998l ) may help to 
constrain the epoch when most primordial disks around high/i ntermediate-niass stars turn 



into debris disks. Cluster memb ership is well determined (e.g. lHerbig||l998l : iLuhman et al. 
19981 . I2OO3I : iMuench et al.ll2007l ). At only 320 pc distant, the cluster provides a sensitive 
probe of low-levels of disk emission in spite of its high infrared background. 



Recently, iLada et al.l (120061 ) and iMuzeroUe et al.l (120061 ) analyzed Spitzer Cycle 1 IRAC 
and MIPS IC 348 data and found evidence for substa ntial disk evolution. Based on the 
IRAC 3.6 fim to 8 fim flux slope, iLada et al.l (120061 ) show that only ^ 30% of cluster 
stars have strong ('thick') IRAC excess emission comparable to stars with primordial disks 
(e.g. most disk-bearing stars in Taurus). Another ^ 23% have weaker ('anemic') IRAC 
excess emission, and ~ 44% lack IRAC excess emission ('diskless'). However, inferring the 
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physical state of disks from the observed IRAC slope is not straightforward. Both debris 
disks and 'evolved' primordial disks ('transitional' disks: those with lower disk masses and/or 
inner holes) can have weak IRAC-excess emission (jCurrie et al.ll2007bl . l2009l ). Furthermore, 
while SED analysis including longer wavelength data (e.g. 24 ^m) helps to break this 
degeneracy (ICurrie et al.ll2009l ). only ?s 30% of IC 348 stars have Cycle 1 MIPS detections. 
Therefore, deeper MIPS data coupled with other diagnostics of the disk evolutionary state 
(e.g. gas accretion signatures) would help to make IC 348 a better laboratory for studying 
disk evolution. 

In this paper, we combine new, high signal-to-noise optical spectra and new, deeper 
MIPS 24 fim and 70 /im photometry of IC 348 with archival data to study disk evolution 
in more detail. We first use the optical spectra to search for evidence of circumstellar gas 
accretion, which is absent in debris disks. Deeper MIPS data is then used to establish a 
much larger sample of disk-bearing stars. From these data, we compare the disks SEDs 
to models to constrain their evolutionary states, to compare disk properties around high- 
mass stars to low-mass stars, and to investigate the timescale for and the morphology of the 
primordial-to-debris disk transition. 



2. New Observations and Archival Data 
2.1. New Optical Spectroscopy of High/Intermediate-Mass IC 348 Stars 



To investigate accretion signatures for IC 348 stars, we took high signal-to-noise spectra 
of cluster stars with published spectral types earlier than KO. Our primary motivation 
for selecting this sample is that these stars lack recent high signal-to-noise observations. 
High signal-to-noise observations are effective in identifying gas accretio n signatures a round 
solar and subsolar-mass cluster stars not previously known as accretors (jDahmll2008l ). New 
observations of early/intermediate type stars may identify sources with emission line reversal 
i n the cores of their Hq, absorption lines, a feature which is consistent with gas accretion 
jDahnJl2008h . 



Selected IC 348 stars were observed on March 1, 2008 with the FA ST spectrograph on th e 
1.5m Tillinghast telescope at Fred Lawrence Whipple Observatory (IFabricant et al.lll998l ). 
The spectra were taken using a 300 g mm~^ grating blazed at 4750 A and a 3" slit. These 
spectra cover 3700-7500 A at 6 A resolution. T he data were processed using the standard 
FAST reduction pipeline (IFabricant et al.lll998l ). The continuum signal to noise at ~ 6550 
A ranged from 80 to 170. To derive line fluxes in H^, we fit the line to a gaussian profile 
with standard IRAF routines. The properties of each source and their spectral features are 
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listed in Table [H 



2.2. Archival Optical/Near-IR Ground-based Data and Archival Spitzer 

3.6-24 fim Data 

We collected publish ed optical and inf rared photometric data and spectroscopic data 
for IC 348 stars from the iLada et al.l (120061 ) catalog. This catalog yields photometric data 
from 3.6 /im through 24 /im, optical extinction (Ay), spectral types, and IR AC SEP slopes 



("thick , anemic o r disk less) . We added BVI photom etry of IC 348 sources from lLuhman et al. 



and bolometric luminosities from Muench et al. 



((20031), iHerbid (119981 ). and iLuhman et al.l (120051) near-IR JHK^ photometry from 2MASS 



(l2007f ). 



Nearly all of the 307 sources in the iLada et al.l (120061 ) sample have spectral types (304) 
and Ay measurements (295). Most have 2MASS data (280), I-band data (270), and IRAC 
8 fim data (265). Completeness in the MIPS 24 band is poorer (88 sources). Many 
sources without MIPS detections are either not detected in the 8 fim IRAC channel (35) 
and/or have photospheric colors through 8 fim (128). 



2.3. New MIPS 24 nm and 70 /xm Photometry 

IC 348 was reobserved in five separate epochs with the Multiband Imaging Photome- 
ter for Spitzer (MIPS) on September 23-27, 2007 (Observing Keys 22961920, 22962176, 
22962432, 22962688, and 22962944; PI J. Muzerolle). The area of coverage is ^ 0.5 square 
degrees. A series of papers focused on analyzing the IC 348 sources with variable mid-IR 
emission will be published later using these MIPS data and IRS spectra (J. Muzerolle, K. 
Flaherty, et al. 2009a,b, in prep.). Our motivation for using these data is to identify new 24 
/im detections of cluster stars from stacking together pointings from the five epochs. Prelimi- 
nary findings indicate that some low-mass IC 348 stars have 24 /im emission varying by up to 
several tenths of magnitudes from epoch to epoch (J. Muzerolle, pvt. comm.). Therefore, by 
combining observations from these epochs together, we sacrifice any information about disk 
variability. The resulting photometry is an epoch-averaged value: photometry for low-mass 
stars from individual epochs may differ by several tenths of magnitudes. 
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2.3.1. MIPS 24 fiTn Image Processing, Photometry, and Source Matching 



Using MOPEX (IMakovoz and Marleaul l2005l ). we first mosaic the Basic Calibration 
Data (BCD data) together using a bicubic sphne interpolation. In constructing the final 
mosaic, we rejected and reinterpolated over pixels with values more than 5 a away from 
the (clipped) mean value. Standard deviation mosaics were also constructed. The Basic 
Calibration Data (BCD) frames have an integration time of 3.67s; most pixels had > 110 
separate pointings. The typical cumulative integration time / pixel i s then ~ 410s, about a 
factor of five greater than the Cycle 1 data analyzed in iLada et al.l (120061 ). 



Source detection and photometry was performed using APEX (IMakovoz and Khan 



2OO5I ) in single-frame mode. In a first pass, we identified sources as clusters of 4-50 pix- 
els in the signal-to-noise mosaic, which is produced by dividing the background-subtracted 
science mosaic by a map of the background fiuctuations (produced from the gaussnoise 
module from APEX). We used the 'peak' algorithm with a 45x45 pixel box to define the 
background and a 2.5cr detection threshold. While these settings yielded robust detections of 
many faint stars, they sometimes 'deblended' single stars, especially those with bright Airy 
rings. Furthermore, these settings missed detecting some sources located in regions with 
smaller-scale fiuctuations in the background (e.g. near the center of the nebula). Therefore, 
in a second pass, we chose a 25x25 pixel box, using a 5cr detection threshold and the 'combo' 
algorithm, which alleviated these problems. 

Pixel Response Function (PRF) photometry was performed on detected sources, using 
an empirical PRF derived from bright (S/N > 50-100) stars in low, uniform background 
regions resampled from the native pixel scale (~ 2.45" /pixels) by a factor of four. The em- 
pirical PRF and model PRF supplied by the Spitzer Science Center show strong agreement. 
PRF fitting was performed over a 5x5 pixel area with fitting tolerance of 10~^. Photomet- 
ric errors were estimated from both the signal-to-noise from background fiuctuations (due 
to both photon noise and nebulosity) and errors in the PRF fitting. We detect ~ 5,900 
candidate point sources; the source counts peak at [24] ~ 11. 

To understand the empirical photometric scatter in our MIPS photometry and to iden- 
tify candidate pre-main sequence stars, we cross-correlated our MIPS-24 fim catalog with 
the 2MASS point source catalog. Using a 4" matching radius, we found 451 matches. The 
left panel of Figure [2] shows the J-K^ vs. Ks-[24] colors for these sources. Overp lotted as a 
dotted l ine is the locus of photospher i c colo rs from the Kurucz-Lejeune stellar atmosphere 
models (jKuruca 1 19931 : iLejeuene et al.l 119971 ) as listed by the STAR-PET tool available on 
the Spitzer Science Center website. Sources with photometric uncertainties in their K5-[24] 
colors < 0.2 magnitudes (grey dots) either lie well redward of the locus (IR excess sources) 
or are tightly clustered along the locus. Sources with larger photometric uncertainties (black 
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dots) appear to have a wider dispersion in color. 

Figure [2] (right panel) shows the [24] vs. Ks-[24] color magnitude diagram of these 
sources. The vertical curves (dotted, dashed, and three-dash dots lines) show the median 1 
and 2 a uncertainties in Kg- [24] colors computed in 1-magnitude bins, assuming an intrinsic 
dispersion of 0.1 magnitudes centered on Ks-[24] = 0.25. There is little blueward deflection 
of colors outside of the 1 a limit; the 2 a limit completely contains all sources with colors 
bluer than Ks-[24] = 0.25. We do not find a significant trend between photometric scatter 
and background level (not shown). Thus, if there is an intrinsic ^0.1 magnitude dispersion 
in colors due to differences in spectral type/reddening, the observed photometric scatter is 
explainable in terms of the computed photometric uncertainties for MIPS and 2MASS. 

Using a matching radius of 4", we identify 27 IC 348 cluster stars with new 24 /xm 
detections. These detections were confirmed by inspection of the final mosaic image, the 
background-subtracted mosaic, and the signal-to-noise mosaic. The coordinates were then 
checked against the positions of other sources to identify possible blends. One source, ID 
9099, is ~ 3" from ID 99. Its catalog position is offset from the MIPS detection by 2.83". 
Another source, ID 1937, is within ~ 4.3" of ID 1928 and has an offset of 3" between its 
catalog position and position on the MIPS-24 /im mosaic. We remove these sources from 
the list, leaving 25 non-blended sources detected using APEX. 

To find mi s sed so urces, we manually searched the positions of the remaining 192 stars in 



the iLada et al.l (120061 ) catalog without APEX detections. This search yielded two additional 
detections (IDs 124 and 373). PSF- fitting photometry was performed on these two sources 
using the IDL version of the DAOPHOT package with empirical PSFs constructed from 
bri ght stars un c ontam inated by nebulosity (S/N > 50-100). With these new data added. 



the iLada et al.l (120061 ) catalog includes 27 sources with new detections, bringing the total 



number of MlPS-detected star s to 115. This new sample is a ~ 30% larger than the sample 



analyzed in iLada et al.l (120061 ). Table [2] lists the properties of the IC 348 stars with new 



MIPS 24 detections. 



2.3.2. MIPS 70 fim Image Processing, Photometry, and Source Matching 

Using MOPEX, we mosaiced both the MIPS-70 fim BCD data and the filtered BCD 
data. As with the 24 /xm data, we combined the individual frames together, rejecting and 
reinterpolating over pixels more than 5a from the (clipped) mean pixel value. We constructed 
the final mosaic images and standard deviation mosaics. The MIPS 70 /im data have an 
intergration time/pixel of 4.19s for the BCD frames and a cumulative integration time of ~ 
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210s, a factor of five improvement over previous data presented in iMuench et al.l (120071 ) 



Source identification and photometry were performed using APEX in single-frame mode. 
For source detection on both the unfiltered and fihered mosaics, we chose a 5a detection 
threshold for initial cluster detection, limiting image segmentation to candidate sources be- 
tween 5 and 250 pixels in area. Because there were too few bright point sources detected 
at 70 fim, we used the model PRF for photometry instead of constructing an empirical one. 
Each source was fit by the model PRF over a 7x7 pixel area with a resampling factor of four 
and a fitting tolerance of 10"^. Photometric errors were calculated in the same manner as 
the 24 fim data. Because sources brighter than ^ 1 Jy can have their fiuxes attenuated on 
the filtered mosaic, we chose the unfiltered mosaic for 70 photometry. 

Because MIPS-70 is not sensitive to stellar photospheric fiuxes of stars at the dis- 
tance of IC 348, we cannot investigate the photometric reliability o f the 70 fim as we did with 
the 24 fim data. We matched the APEX 70 /im source list to the iLada et al.l (120061 ) catalog 
using a 5" r natch ing radius. We find 6 sources in the MIPS-70 catalog with matches in the 
Lada et al.l (120061 ) catalog, four of which (IDs 31, 67, 10343, and 10352) are new detections. 
We visually confirmed these detections on both the unfiltered and filtered mosaics. Table [3] 
lists the properties of stars detected at MIPS 70 fim. 



2.3.3. MIPS Upper Limits for Unmatched Sources 



Figure [3] shows a RGB mosaic image produced fr om our 24 iirn and 70 //m mosaics and 
the post-BCD mosaic of the IRA C 8 nm data used by lLada et al.l (120061 ). Most of the IC 348 
sources in the iLada et al.l (120061 ) catalog are located within regio ns of enhanced ne bulosity. 
While nebulosity does not serious impede detections of most of the lLada et al.l (120061 ) sources 
at 8 /im, the decreased sensitivity of MIPS and enhanced nebular emission at 24 //m and 70 
/im greatly increase the number of nondetections at these wavelengths. 

To compute the 5o" upper limits for sources lacking detections, we use the 24 /xm and 70 
/im uncertainty images as maps of the noise. T he median ba fiux density is computed over 
a 3x3 region centered on the coordinates of each iLada et al.l (120061 ) source. To calculate the 
5a" fiux upper limit for MIPS 24 /xm, we multiply the fiux density by a pixel area of 7r(2.45 
pixels)^ and use an aperture correction of 1.697. For the 70 /im data, we multiply the fiux 
density by a pixel area of 7r(1.62 pixels)^ and use an aperture correction of 1.735. Upper limits 
at 24 /xm varied from k, 0.3 mJy ([24] ~ 11.0) in low (black) background regions to ~ 180 
mJy ([24] ~ 4) in the high background regions of the nebula center (reddish- white/white). 
Typical upper limits at 70 /xm varied from ^ 30 mJy in low background regions to ~ 420 
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mJy in the center of the nebula. 



2.3.4- Comparisons with Previous IC 348 MIPS Photometry 



In Figure HI we compare our MIPS photometry with results from lLada et al.l ( 20061 ) and 
Muench et al.l (120071 ) derived from Cycle 1 data. The left panel contrasts our [24] magnitudes 
and upper limits, where grey circles represent magnitudes for IC 348 members detected from 
both data sets and crosses represent 2a upper limi t s for I C 348 members lacking detections 
from both data sets. Compared to the iLada et al.l (120061 ) photometry, our [24] photometry 
exhibits a systematic offset of -0.0075 magnitudes and a median dif ference of 0.05 niagnit udes. 
While our 70 flux for ID- 13 is just over half that quoted bv iMuench et al.l (120071 ) . our 



measurement for ID-51 is nearly identical to that from lMuench et al.l (120071 ). 



In the cases where the [24] photometry differs by more than 0.2-0.3 magnitudes, there is 
evidence that our photom etry is more accur ate. For example, ID-7 has [24] = 7.61 from this 
work and [24] = 6.95 from lLada et al.l (120061 ). Comparing its flux densities on our mosaic and 
the Cycle 1 data post-BCD m osaic indica t es tha t variability is unlikely to explain the strong 
disagreement in photometry. iLada et al.l (120061 ) find that ID-7 is about equal in brightness 
to nearby ID-8. However, it is clear that ID-7 is much fainter because its peak flux density 
above the local background is just over half that of ID-8's (~ 3.7 MJy sr~^ vs. 7.1 MJy sr^^). 
Thus, ID-7 should be ~ 0.7 magnitudes fainter than ID-8, as confirmed by our measurements 
([24]/z)-7=7.61 and [24],^_8=6.89). 

ID-19 is another source whose dis crepant photomet ry probably cannot be due to vari- 
ability ([24]=6.24 here, [24]=4.6 from iLada et al.l l2006l ). This source is located near the 



center of the nebula, where the background flux density increases from ^ 75 MJy sr"^ to ^ 
98 MJy sr~^ over ^ 20". Unless we assume that the average background is ~ 75 MJy sr~^, 
we fail to reproduce the iLada et al.l (120061 ) result from either PRF fitting or PSF fitting. 
However, there is a significant slope in the local background at constant declination near 
ID-19; thus, the average or median background is much higher than 75 Myr sr~^. Indepen- 
dent MIPS 24 fim photometry yields values for ID-7 and ID-19 that are similar to ours (J. 
Hernandez 2009, private communication). 

The left panel of Figure HI also compares our upper limits for nondetections (crosses). 
Assuming the iLada et al.l (120061 ) upper limits are accurate. The dotted line corresponds to 
the predicted locus of upper limits given our integration time in the photon noise limit 
(2.51ogio(v^) ~ 0.87). Because nebular emission dominates the background in most regions, 
our upper limits should lie to the left of the line. However, most of our upper limits are 
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clearly fainter than predictions based on Lada et al.'s upper limits. Independent analyses 
also find significantly fainter upper limits (J. Hernan dez 2009, private communication). The 
disagreement likely arises because iLada et al.l (120061 ) add together the background flux and 
its uncertainty, which would determine the upper limit of the source plus background, not 
the upper limit of the source flux, since the background is an astronomically distinct object. 

Even though the new MIPS data are deeper, the ri ght panel o f Figur e HI shows that our 
photometric uncertainties are larger than those listed by lLada et al.l (120061 ). Our results agree 
with those of Lada et al. much better if we assume that the background uncertainty (from 
nebu lousity) doe s not c ontribute to photometric uncertainty. In spite of some disagreements 
with iLada et al.l (120061 ) 's photometry, results of the analysis presented in the proceeding 
sections are not qualitatively affected by which photometric catalog is used for previously 
detected stars. 



Results from New Data 



We now use new data to expand on Lada et al's investigations of the IC 348 disk 
population. First, we investigate accretion signatures to probe the presence or absence of 
warm circumstellar gas in disks surrounding cluster stars. We then analyze the mid-IR 
colors of stars with ne w MIPS dete c tions to identify IR excess sources. We then follow 
the analysis methods of ICurrie et al.l (120091) to model the SE Ds of new sources to constrain 
their evolutionary states and to refine the iLada et al.l (l2006l ) analysis of previously detected 



sources. 



3.1. Gas Accretion Signatures and Upper Limits 

To investigate the gas content of high/intermediate- mass stars, we consider line fluxes 
derived from new, higher signal-to-noise optical spectra. Our spectra show a variety of Hq, 
morphologies (Figure [T|). IDs 3, 8, and 30 lack 2a fluctuations that could be interpreted 
as Hq, emission in the central absorption core. The spectra of ID-2 shows a flattening of 
the Hq absorption core, though it lacks unambiguous evidence for emission line reversal. 
Higher-resolution spectroscopy s uggests that t his source may have emission line reversal (S. 



Dahm, unpublished) as noted in iDahml (120081 ) 



ID-6 has Ho in absorption but a smaller measured equivalent width than the typical 
width for early/mid G stars (~ 3-4 A). Thus, the H„ line may be filled in by emission. 
However, the spectrum shows no clear evidence for an emission line core. Furthermore, 
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ID-6 is very chromospherically active: it has the larg est X-ray luminosity of any Chandra - 
detected source in IC 348, ~ 1.96x10^^ ergs s~^ (jPreibisch and Zinneckerl l200ll . |2002| ). 
Such activity should produce weak Hq, emission and result in a weakened absorption line. 

ID-21 shows Ha in emission with EW(Hq,) ~ 5A consistent with previous results. Sur- 
prisingly, ID-19, an A2 star not previously known to be accreting, also shows evidence for Hq, 
emission consistent with accretion. Even though the line is in absorption, a central emission 
core is clearly visible. 

For the sources lacking clear evidence for accretion, we can rule out Hq, line emission 
hiding in the nois e at a l evel o bserved in late- type stars. We estimate the upper limit of the 
Hq, flux following lAUenl (119951 ). First, from the spectral type and extinction, we derive the 
expected dereddened continuum flux in the Cousins R band. We then measure the equivalent 
width of the continuum noise near the Hq line. An upper limit to the line flux then results: 
FluxH^^ui ~ EW(noise)xR/;„^. 

Finally, we derive an upper limit to the accretion rate scaled to that of I D-21 from Lff. 



( X GM M^/R^. We determine R^, from the de ri ved v alues for Tg and from lMuench et al. 



(120071 ) and use the iD'Antona and Mazzitellil (119941 ) isochrones t o estimate stellar masses 



Although the theoretical ages of young s tars are uncertain (e.g., ID'Antona and MazziteUi 
I994J : ISiess et al.ll2000l : iBaraffe et al.lll998l ) the masses derived from and Tg using different 
isochrones are generally consistent. As shown in Table 1, the Hq, line flux (accretion rate) 
for sources lacking evidence for gas accretion is < 10% (< 0.5-10 %) of the Hq emission 
(accretion rate) in ID-21. These upper limits also rule out accretion at a rate equal to that 
for ID-19. 

As with any pre-main sequence star, it is still possible that the stars lacking clear 
accretion signatures are accreting at undetectable rates. Furthermore, it is possible that 
disks surrounding these stars lack warm circumstellar gas accreting onto the star but contain 
some gas in outer disk regions. Their optical spectra and resulting upper limits on accretion 
are more consistent with a lack of accretion and thus a heavily depleted or even nonexistent 
gas reservoir. 



3.2. Mid-to-Far IR Colors of IC 348 Stars with new MIPS Detections 



3.2.1. MIPS 24 fxm Colors 



Figure O shows the [3.6]-[8] vs. [8]-[24] color-color diagram for new MIPS detections. 
Overplotted are the la error bars in each filter. We include a reddening vector of length 
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Ay =20 (the maxi mum extinction), using the extinction laws from llndebetouw et al.l (120051 ) 
and iMathid (Il990l ). Most newly-detected stars have spectral types later than KO (left panel; 
grey dots); only one new detection has an earlier spectral type (left panel; black dot). 

Cluster stars with new 24 fim detections generally have excess emission consistent with 
reprocessed radiation from circumstellar disks. Most new detections have an observed [3.6]- 
[24] color redder than ~ 1.5-2 magnitudes. In contrast, [3.6]-[24] < 0.5 for bare s tellar 
photospheres for stars of most spectral types (e.g. Currie. Plavchan. and Kenyon 20081 ). Be- 
cause the maximum extinction for cluster members of Ay ~ 20 (shown as an arrow) reddens 
the stars by less than ~ 0.5 magnitudes, high extinction is unlikely to explain these colors. 
Therefore, the excess emission plausibly arises from circumstellar material. 

The distribu tion of mid- I R colo rs is consistent with results for previously detected cluster 
stars an alvzed inlLada et al.l (120061 ). Most stars identified as 'diskless' (right panel; light grey 



dots) by Ihada et all J2006h either have photospheric MIPS emission or excess emission less 



than 2 magnitudes. Their MIPS excess emission is less than that from stars with strong 
('thick') IRAC excess (right panel; black dots) and stars with weak ('anemic') excess (right 
panel; dark grey dots). Both 'anemic' and 'thick' disks have a range of [8]- [24] colors from 
~ 2 to 4.5 and differ primarily by their levels of IRAC emission. 

To compare the mid-IR colors of IC 348 sources with fiducial values, we overplot colors 
of the median Taurus SED (asterisk) reddened to E(B-V) ~ 1.56, the mean v alue for IC 348 
stars. To produce the Taurus colors, we use fluxes fro m iFurlan et al. I J2006h for the 8 fim 
and 24 fim filters and approximate the 3.6 /im flux from lD'Alessio et al.l (119991 ). The Taurus 
colors lie within the distribution of colors for 'thick' sources in IC 348. 

While many newly-detected stars have levels of excess emission typical of primordial 
disks, the full range of IRAC vs. MIPS colors hints at a diversity in disk properties. Two 
stars have nearly zero [3.6]- [8] color and [8]- [24] > 3. These colors are consistent with stars 
whose disks have inner holes cleared of dust and are optically-thick at larger stellocentric 
distances. A second group of stars also have zero IRAC color but may have weak 24 fim 
excess emission ([8]-[24] ~ 0.5-2). These colors are consistent with the colors of either debris 
disks or remnant protoplanetary disks that are actively disappearing. 



3.2.2. MIPS 70 fim Colors 

Figure [6] shows the distribution of [24]-[70] colors for stars detected at 70 fim. All of 
these stars have spectral types later than KO and have substantial emission at both 24 fim 
and 70 fim. The right panel of Figure M shows that the stars with 70 detections include 
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both 'anemic' and 'thick' sources, indicating that at least some stars with dust-poor inner 
regions may have very luminous outer disks. While the distribution of MIPS colors appear 
correlated (stronger 24 /im emission — stronger 70 fim emission), the small number of stars 
with 70 fim detections precludes any robust statistical analysis to determine if this trend is 
real. 

To provide a context for our 70 fim detected sta rs, we compare their colors to those for 



two stars in the 77 Cha Association using data from iGautier et al.l (120081 ). Even though r] 
Cha is substantially older than IC 348 (~ 8 Myr), its proximity to the Sun and lack of high 
IR background allow for many 70 fim detections; the median Taurus SED currently does not 
include a 70 fim data point. Overplotted as diamonds are colors for ECH-15 (top center), 
an accreting late- type star with a primordial disk, and ECH-2 (lower left), which has weak, 
optically-thin emission consistent with a debris disk. The MIPS colors of IC 348 stars are 
slightly redder than ECH-15 at both 24 fim and 70 fim and far redder than ECH-2. 



3.3. SEDs of new MIPS Detections 



3.3.1. Constraints on Disk Properties of Sources with MIPS-70 fim Detections 



The addition of the MIPS-70 fim data better determines plausible circumstellar disk 
properties from models with nearly identical fluxes through 24 fim but very different far- 
IR fluxes. To better constrain the disk morphology of IC 348 sources, we compare the 
SEDs of MIPS- 7 /xm-detec ted sources with the grid of 3D radiative transfer models from 
Robitaille et al.l (120061 . 120071 ). In fitting the SEDs, we assume a 10% uncertainty in the IRAC 
and MIPS-24 fim fluxes and a 20% uncertainty in the MIPS-70 fim fluxes to account for 
photometric error, zero-point uncertainty, and variability. We allow for a range of distances 
(~ 300-350 pc) and a range of optical exti nction, typical l y cora parable to the standard 
deviation in the extinction as determined by iMuench et al.l (120071 ). The online fitting tool 
then returns the SEDs of the best-matching models and lists the disk characteristics (inner 
and outer disk radius, accretion rate, disk mass, disk flaring index, etc.) of the models. 



Figure [7] shows the fitted spectral energy distributions. In each panel, we overplot the 
SED of the best-fit bare stellar photosphere (dotted line), the best-fit SED model (solid black 
line), and other SED models with x^-xLst < ^O- The xLs* sources range from x^=9.4 

(ID 103542) to = 136 for ID 67. While the value for the latter source is large, much of 
the power in the values comes from measurements at shorter wavelengths (e.g. J band) 
where the photometric uncertainty is very small. 



One source, ID-51, has a nearly flat SED (AFA) from 2 fim through 70 fim and is thus 
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consi stent with "flat-spectrum" protostars like those found in Taurus- Auriga (IKenvon and Hartmann 



1995 ). Such sources may be in the process of dispersing their protostellar envelopes (ICalvet et al 
19941 ) ■ Parameters for the best-fit model include a flaring index of (3 =1.1 {(3 = H/r, where 
H is the disk scale height) and a disk mass of ~ 3 x 10~^ Mq. The SED of ID-51 rises 



from I band through 5.8 ^m, which is consistent with a 360 K protostar as suggested by 
Muench et all J2007h . 



Half of the 70 sample - IDs 13, 10343, and 10352 - have SEDs consistent with pre- 
main sequence stars surrounded by primordial disks. These sources have clear IR excesses 
at near-IR wavelengths and SEDs whose flux declines with wavelength from peak values 
at ~ 1-2 fim. Assuming a standard grain size distribution (n(a) oc a~^'^) and a solar gas- 
to-dust ratio, these sources have disks with masses ~ 10^^-10^^ Mq, comparable to the 
Minimum Mass Solar Nebula. The disk flar ing indices range from /3 = 1.01 (for ID 13) to 
1.17 (ID 10343). While iMuench et all J2007h identify ID-13 as a protostar, its IRAC slope is 
slightly steeper than their criteria for either protostars or flat spectrum sources. Our analysis 
coupled with its published spectral type (MO. 5; Tg ~ 3800 K) indicates that it may be at the 
boundary between a flat spectrum source and a pre-main sequence star + primordial disk. 

Two sources, IDs 31 and 67, show evidence for a strong flux deficit from 3.6 fim through 8 
fim but substantial emission at both 24 fim and 70 fim. Disk parameters from the best-fitting 
models have a wide range of disk masses (10~^ - 10~^ M©) but typically have large inner disk 
boundaries (~ 1-50 AU). While a depleted mass in warm dust and inner holes are consistent 
with disks in an intermediate stage between primordial and debris -'evolved' primordial disks 
or, as typically described, 'transitional' disks- there are alternate possibilities. Specifically, 
given that IC 348 contains many protostars, the strong 24 and 70 /zm emission could be due 
to a protostar lying in close proximity to IDs 31 and 67. Because such sources would be 
unresolved with MIPS and too faint to detect from optical/near-to-mid IR imaging, their 
presence cannot be ruled out. However, t hese two sourc e s are located in regions well away 
from most IC 348 protostars identified by lMuench et al.l (120071 ). 



3.3.2. Evolutionary States of Disks with New MIPS-24 /iw, Detections 



In this section, we investigate the evolutionary states of disks surrounding sources with 
new MIPS-24 fim detections. To provide a simple but effective probe of the disks' evolution- 
ary states, we compare their SEDs to the median Taurus SED and upper/lower quartiles of 
Taur us sources and debr is disk SEDs. We use the Taurus SED and upper /lower quartiles 
from iFurlan et al.l (120061 ): for the debris disk SED, we adopt the synt hetic SED of a warm, 
terrestrial planet-forming debris disk around a solar-mass star from iKenyon and Bromley 
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( 120041 ) ■ Evolutionary states derived from these simple cor nparisons show go od agreement 
with those obtained from more detailed SED modeling (see lCurrie et al.ll2009l ). 



Following ICurrie et al.l (120091 ) . we classify sources as having primordial disks if their 3.6 
fim to 24 fim emission is comparable to or greater than the lower quartile of the median Tau- 
rus SED. Sources lacking clear excess emission in the IRAC bands but emission comparable 
to Taurus at 24 fim are identified as disks with inner holes. Sources with IRAC and MIPS 
excess emission significantly weaker than the lower quartile T aurus SED but str onger than 
the debris disk SED are labeled 'homologously depleted disks' (ICurrie et al.ll2009l ). Here, we 
define "significantly weaker" as less than half of the lower -quartile Tau r us SE D flux at 24 
microns. This stricter criterion is justified because, unlike in lCurrie et al.l (120091 ) . most (68%) 
of the new detections are stars later than M3, which are intrinsically redder and fainter than 
the K7-M2 stars used to derive the median Taurus SED, and thus should have systematically 
weaker primordial disk emission. Finally, sources with IR excess emission comparable to the 
debris disk SED are identified as debris disk candidates. 

Figure M shows SEDs of sources representative of each disk class. All source SEDs 
are labeled by their ID number, J2000 coordinates, and extinction. We also include the 
IRAC slope (st ron/weak IRAC, ph otosphere) following the formalism of iLada et al.l (120061 ) 
as modified by ICurrie et al.l (120091 ). The figures are arranged to illustrate the evolution of 
primordial disks (top panel) into debris disks (bottom panel) along two separate pathways: 
homologously depleted disks (second and third row panels, lefthand side) and disks with 



inner holes (second and third row panels, righthand side) (see lLada et al.ll2006l : ICurrie et al. 



20091 ). The Appendix contains the atlas of all SEDs; Table [5] lists results for all stars. Below 
we describe the frequency and characteristics of each disk class. Because of a lack of short 
wavelength data, we did not determine the evolutionary states for IDs 746 and 2096. 



1. Primordial Disks- Nine sources have SEDs indicative of primordial disks. The levels 
of disk excess from most of these sources are more comparable to the lower quartile of 
the median Taurus SED than the median Taurus SED. However, most sources in IC 348 
have spectral types later than M4; it is not clear whether their sli ghtly weaker eraission 
is an evolutionary effect or simply a temperature (contrast) effect (lErcolano et al.ll2009l ) 
because the median Taurus SED was derived from hotter (K7-M2) stars. Some of these 
sources have mid-IR emission comparable to the median Taurus SED, a characteristic 
consistent with a relatively unevolved primordial disk. 



2. Homologously Depleted Disks- Seven sources have disks that produce a clear IR 
excess from ~ 4.5 /im to 24 /xm but at a level significantly weaker than Taurus: a 
'homologously depleted' disk (ICurrie et al.ll2009l ). Because both the IRAC and MIPS 



- 15 - 



emission is depleted by up to a factor of 10 relative to the Taurus SEDs, it is unlikely 
that this weaker emission is simply due to flatter disk geometries. The weaker emission 
is likely due in part to signiflcantly lower warm dust masses which can ernerge f rom a 



combination of grain growth and accretion onto the star (see lCurrie et al.ll2009l ). The 



SEDs exhibit a power-law decline with AFa (x X where q ~ 1.5-2.5. 

Disks with Inner Holes - Two sources show evidence for a disk with little or no 
emission in the IRAC bands but substantial MIPS emission, up to a level comparable 
to the median Taurus SED. These properties are characteristic of disks with inner 
regions mostly cleared of material. Because these disks represent only one of the two 
pathways from pri mordial disks to debris disks, and thu s comprise one of two types of 



transitional disks, (ICurrie et al.ll2009l : iLada et al.ll2006l ). we do not equate them with 



transitional disks but call them 'disks with inner holes' instead. 

4. Debris Disks Candidates - Five sources have MIPS emission comparable to pre- 
dicted fluxes from debris disks. This emission is a factor of ~ 100 weaker than pri- 
mordial disks and is also weaker than emission from homologously depleted disks, and 
disk with inner holes. We label these sources 'debris disk candidates'. 

Debris disks are usually deflned as disks with second-generation dust, a property that 
correlates with broadband SEDs but cannot be directly inferred from SEDs. Disks 
around late type, low-mass stars have a much lower contrast with the photosphere 
than do disks around high-mass stars; thus, there is less difference in broadband fluxes 
between debris disks and disks in earlier evolutionary stages. Some of these disks 
could be very evolved (homologously depleted, inner hole) disks with flrst generation 
dust, not second generation dust. Others have larger photometric errors at 24 fim and 
could be bare photospheres instead. Thus, our labeling of these sources as 'debris disk 
candidates' is tentative and requires more analysis for conflrmation. We consider these 
sources further in §3.4. 

5. Stars Without Disks — Two sources have SEDs that are consistent with bare stellar 
photospheres. ID-47, a K4 star, is the source that most obviously lacks IR excess 
emission. It lacks any IRAC-excess emission, and its 24 /xm flux is slightly bluer than 
a stellar photosphere. However, this source has a marginal, 5 a detection, making 
its flux consistent with a stellar photosphere to within ~ 2 a. ID-49 also has a very 
marginal IR excess at 24 /im (Ks-[24] ~ 0.5) that is within 2 a of photospheric colors. 
Thus, there is no clear evidence for a disk surrounding this star. 
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3.3.3. Mid-IR Colors of Sources with Disks in Various Evolutionary States 

In Figure [9] (top panel), we show the dereddened Ks-[8] vs. Ks-[24] color-color diagram 
for sources with new MIPS-24 /im detection^. Symbols for sources with new MIPS de- 
tections correspond to their derived disk evolutionary state: primordial disks (circles), disks 
with inner holes (triangles), homologously depleted disks (inverted triangles), and debris disk 
candidates/photospheres (squares). Sources with primordial disks lie redward of a region de- 
fined by a vertical and horizontal short-dashed line. Polygonal regions enclose the range of 
colors for the two types of transitional disks: homologously depleted disks (dash-three dots) 
and disks with inner holes (dashes). Debris disk candidates and photospheres comprise 
all the newly-detected sources bluer in Ks-[24] than the homologously depleted/inner hole 
boundaries. The positions of sources later than KO in IRAC/MIPS color-color diagrams thus 
show a good separation by disk evolutionary state as determined in the previous section. 

For reference, we overplot the colors of median Taurus SED with upper and lower 
quartiles (asterisks). The median Taurus SED (Ks-[8]=2.1, K5-[24]=5.1), the upper quartile 
(K,-[8]=2.5, K,-[24]=6.3), and lower quartile (K,-[8]=1.8, K,-[24]=4.8) all lie within the 
primordial disk region. All of the sources identified as having primordial disks have 8 (24) 
/im fluxes (measured relative to K-band) a factor of > 1.6 (2.1) smaller than the lower 
quartile Taurus SED or 24 fluxes reduced by a factor of ~ 2.1. The population from 
which the median Taurus SED was constructed includes sources with inner holes (e.g. DM 
Tau) or a significantly low dust mass (V836 Tau), so a median SED constructed only from 
sources with full disks has redder colors. Therefore, the fact that the boundaries enclosing 
primordial and non-primordial disks in Figure [9] (top panel) are bluer than positions of the 
Taurus SED colors indicates that our identification of non-primordial disks is conservative. 



We also overplot IC 348 sources with previous MIPS detections from lLada et al.l (120061 ) 
as filled dots with a brightness corresponding to their IRAC slope: 'thick' sources (light-grey 
dots), 'anemic' sources (grey dots), and 'diskless' sources (black dots). There is generally 
a good correspondance between the IRAC slope and the inferred evolutionary states as de- 
fined by the colors (e.g. diskless=photospheres/debris, anemic=homologously depleted/inner 
holes, thick=primordial) . However, some 'thick' sources are located in the region correspond- 
ing to homologously depleted disks and some 'anemic' sources are located in the region 
corresponding to primordial disks. 

Figure [HI (bottom panel) shows the dereddened Ks-[8] vs. Ks-[24] diagram for sources 
earlier than KO, which for a 2-3 Myr-old cluster corresponds to stars with masses > 1-2 



^Two sources - IDs 746 and 2096 - lack 2MASS detections and are thus not shown. 
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Mq flBaraffe et al.lll998l : iD'Antona and Mazzitellilll994f ). There is no fiducial SED (and thus 
no empirical colors) characterizing primordial disks surrounding stars of this mass range. 
To provide a context for the colors of these IC 348 sources, we ove rplot the colors of a 
F9 star with a massive debris disk, h and x Per-S5 (small asterisk, ICurrie et al.l l2007bl ). 
and AB Aurigae (large asterisk), an AO star in Taurus surrounded by a primordial disk. 
The colors for A B Aur are estimate d from the 2MASS and Spitzer-IRS fluxes at 2.2 fim, 8 
/im, and 24 fim (IFurlan et al.ll2006l ). ID-4 has colors much bluer than h and x Per-S5 and 
is thus consistent with being a debris disk candidate. Among previously-detected sources, 
one (ID-8012) clearly has primordial disk emission (Ks-[8] ~ 3.4, Ks-[24] ~ 5.3). Five others 
with Ks-[24]=2.5-5.5 require further analysis because their mid-IR colors could be consistent 
with either a debris disk or a remnant protoplanetary disk (i.e. homologously depleted, inner 
hole) . 

In summary, the dereddened Ks-[8] vs. Ks-[24] color-color diagram separates late-type 
IC 348 sources by disk evolutionary state. This analysis is more ambiguous for early /intermediate 
spectral type stars. Sources of any spectral type identified as debris disk candidates on the 
basis of their mid-IR colors require further analysis because debris disks, by definition, lack 
gas accretion and have second generation dust. These properties cannot be determined from 
mid-IR colors alone. In the next section, we will investigate the properties of sources with 
ambiguous states, including all debris disk candidates. After determining their likely evolu- 
tionary states, we will use the results from color-color diagrams presented in this section to 
analyze the global properties of the IC 348 disk population. 



3.4. Evolutionary State of Disks Surrounding Early /Intermediate- Type Stars 

and All Debris Disk Candidates 

Many early /intermediate spectral type sources with mid-IR colors bluer than those 
harboring primordial disks (e.g. AB Aurigae) have ambiguous disk evolutionary states. 
Many sources surrounding late type, low-mass stars with new MIPS-24 detections have 
extremely weak mid-IR disk emission that could either come from a debris disk or a remnant 
protoplanetary disk. In this section, we constrain the evolutionary state of disks surrounding 
stars from these two populations. We first use new and archival optical spectra to identify 
accretors and to remove them from the population of debris disk candidates. For remaining 
debris disk candidates, we estimate the mass of emitting dust (through 24 /im) and calculate 
dust removal timescales. The early /intermediate type stars investigated in this section are 
IDs 2-8, 11, 19-20, 25, and 30-31. The late-type stars investigated here, all of which have 
new MIPS-24 fim detections, are IDs 39, 53, 56, and 1939. 
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3.4- 1- Accretors 

Using new and archival optical spectra, we can remove several accreting early-type stars 
from the list of debris disk candidates. These include ID 19 - identified as an accretor from 
this work - and IDs 5 and 31, which were previously identified as accretors. None of the 
late-type stars show unambiguous evidence for accretion based on their Hq, equivalent widths. 

One source -ID 2 - is a borderline case. Our spectra for ID-2 may exhibit slightly 
weakened Hq, absorption line compared to non-accreting stars of the sanie spe ctral type (ID- 



8). Based on unpublished high- resolution echelle spectroscopy, iDahml (120081 ) identifies this 
star as having an H^, emission line reversal plausibly due to accretion. Thus, we remove 
it from our sample of debris disk candidates. We treat this star as lacking evidence for 
circumstellar gas, though we caution that accretion may persist at low rates. Given that the 
presence or absence of accretion is a primary marker of disk evolutionary state, the inability 
to rule out or identify accretion for this source renders its disk status uncertain. 



3.4-2. SED Modeling of Debris Disk Candidates: Estimating Dust Masses and 

Temperatures 



Next, we compare the SEDs of early-to-late type debris disk candidates with 3D radiative 
transfer disk models for a range of disk masses. To model the near-to-mid IR emission of 
debri s disk candidates , we co mpare their SEDs to the grid of 200,000 radiative transfer models 
from iRobitaille et al.l (120061 ) . which include a range of spectral types and dust masses. The 
infrared fluxes predicted for disks are sensitive to the dust masses down to ~ 10~^° Mq. By 
selecting the best-fitting SEDs, we estimate the minimum mass of small (< 100 /im - 1mm) 
dust grains required to reproduce the observed disk emission. We compare this estimate to 
typical dust masses in debris disks and pri mordial disks. Typical dust masses for debris disks 
range from 10"^^ M© to 5x10"^ Mq (cf. Ichen et allboosl : ILow et al.l[2005[ l. Dust massefl 
typically needed to reprod uce the IRAC and MIPS excess emission for transitional disks 
range from lO'^-lO^^ Mq dCurrie et allboooh . 



We choose a representative sample of sources to model: IDs 6, 8, 20, and 30 for the 
early /intermediate- type population and IDs 39 and 1939 for the late- type population. Be- 
cause most of the other sources have weaker disk emission or comparable emission, the 



^This mass is not the same as the disk masses , which is more typicahy cited. To derive the disk mass 
from the dust mass, Andrews and Wilham"^ ( 2005 ) and other other authors assume a gas to dust ratio of 
100. The Robitaille et ahl ( 2006h models assume the same gas to dust ratio. 
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required mass of emitting dust in their disks is most likely comparable or lower. Therefore, 
modeling the SEDs of these sources to estimate the mass of small, emitting dust grains is 
sufficient to yield approximate dust masses for all debris disk candidates. 

Figure [10] shows SED model fits t o the six sources. Ta ble [6] summarizes the model 
parameters for each source. The best-fit iRobitaille et al.l (120061 ) models accurately reproduce 
the observed fiuxes fr om these sources. Assuming a gas-to-dust ratio of 100, the required disk 
masses for the best-fit IRobitaille et al.l models are ~ 2-8 x 10"^ Mq. The extremely low disk 
masses are consistent with the disks being very optically thin. The required masses in dust 
are ~ 10~^*'-10~^ M©. These dust masses are lower than typical dust masses of primordial 
disks but similar to debris disk masses. Unless these sources had enormous reservoirs of 
extremely cold dust, which would not significantly contribute to the observed IRAC and 
MIPS 24 fim emission, the total mass of dust in small grains is also likely much less than 
typical dust masses for primordial disks. 

We can place further constraints on the temperature and luminosity of the dust. Most 
debris disk c andidates have clear ex cess emission only at 24 /im with dust temperatures colder 
than 200 K (jCarpenter et al.ll2009[ ). However, IDs 6 and 8 have 8 fim excess emission, which 
indicates that their dust, if confined to a narrow ring like in a debris disk, is warmer. Indeed, 
a single temperature blackbody fit to the SEDs of these sources yields dust temperatures 
of 400 K and 270 K. These temperatures are similar to those for know n terrestrial planet- 



forming debris disks such as h and y Per-S5, EF Cha, and HD 113766A (jCurrie et al. 



2007b 



Rhee et al.l l2007al : iLisse et al.ll2008l ). The fractional infrared luminosities (L^j/L^) of debris 
disk candidates from one and two-temperature blackbody fits range from 5x10"^ (ID-6) to 
10~^ (ID-3). These fractional luminosities are two to three orders of magnitudes less than 
those for primordial disk-bearing stars but very similar to the stars wi th warm, terrestria l 
planet-forming disks as well as those with colder debris disks (e.g., (3 Pic lRebull et al.ll2008l ). 



3.4-3. Dust Dynamics and Removal Timescales 



To provide better constraints on the evolutionary state of debris disk candidates, we 
invest igate dust dynamics. If the dust must be second generation, then the disks are debris 
disks (IBackman and Parescdl 19931 ). We determine grain removal timescales in optically thin 
conditions in two cases: where there is a significant amount of undetected circumstellar gas 
(1-10 M©) and where there is little or no gas. If the timescale for removal is much less than 
the age of the star in either case, the grains require a replenishment source through collisions 
between larger bo dies. Disks around such star s would then fulfill the standard definition of 
a debris disk (e.g. 



Backman and Parescdll993l ) 
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Case 1: Optically-Thin Disk with Residual Circumstellar Gas 

Takeuchi and Artymowicz J200lh extensively investigated the dynamics of circumstel- 
lar dust in optically-thin disks with residual gas. Under these conditions, the dust 
grains migrate outward at a rate primarily determined by 13, the ratio of the force from 
radiation pressure to the gravitational force: 



3L^ < Qprja) > 
IGnGM^csps 



(1) 



When (3 > 0.5 the dust is removed from the system entirely on orbital timescales, 
which are ~ 100-1000 yr (see also IStrubbe and Chiangl 120061). Removal occurs in - 



depen dent of whether there is residual gas left in the disk (ITakeuchi and Artymowicz 
200ll . Takeuchi 2008, pvt. comm.). Equation [1] can be rearranged to find the size of 
dust grain below which radiation pressure rapidly removes du st from the system (see 
also iBackman and Parescd Il993l : iBurns. Lamy and Soteii Il979l ) : 



^max,t,m < 1.14/im( — )( 



Igcm 

Ps 



-3 



(2) 



Typically, the sizes o f mid-IR emitting d ust grains range from f» 1 iim to 10 //m (e.g. 
Augereau et al.lll999l : iFurlan et al.ll2006l ). 

For early-type stars listed in Table O Smax ranges from 2.2 to 17.8 pm, assuming 



< Qpr > ~ 1 and p. 



g cm 



Using a more realistic dust grain density of 



0.1 g cm ^ (e.g. lAugereau et al.lll999l ) essentially removes all dust grain sizes that can 



plausibly produce IRAC and MIPS disk emission. Grains responsible for producing 
this emission are blown out; mid-IR emission ceases unless the new grains are produced 
from collisions between larger bodies. 

Dust grains much larger than s^ax {P < 0.5) can produce 8-24 pm broadband excess 
emission in these early type stars and must also be second generation. For a 1-10 
gas disk around an early-type star, the combined influence of radiation pressure, 
Poynting- Robert son drag, and g as drag push 10-100 urn-sized dust from 10 AU to 
~ 300 AU in ^ 600-60,000 yr (ITakeuchi and Artymowica l200ll ). For the early-type 
stars in Table El 120-200 K dust is located at ^ 5-35 AU. If pushed to ~ 300 AU, 
the dust temperature drops by up to a factor of 5, making the dust too cold to emit 
strongly at 24 pm. If new dust at 5-35 AU is not produced, then the disks' 8-24 pm 
emission will become undetectable. Thus, the observed mid-IR emitting dust must be 
second generation. 

The effect of radiation pressure on dust surrounding late-type stars is less clear. The 
maximum grain size for ID-39 is ~ 1.5 pm, a value that is slightly smaller than Smax for 
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the earlier-type stars. However, Smax for ID- 1939 (and, by inference, all M stars that are 
debris disk candidates) is ~ 0.98 fim, which begins to approach the sizes of large grains 
in the interstellar medium. While radiation pressure may be able to remove grains with 
these sizes, it is probably ineffective at completely removing grains with sizes typical 
of circumstellar disks (up to 10 /xm) unless the grains are extremely porous. Therefore, 
the dust responsible for 24 /zm disk emission is not necessarily removed from the disk 
by radiation pressure. However, grains with sizes typical of circumstellar disks and 
and densities of ~ 1 g cm~^ have up to 0.3 and could be transpo rted from their 



original locations to ^ 300 AU scales (ITakeuchi and Artymowical200ll ). 



Case 2: No Gas or Trace Amounts of Gas 

In the absence of residual circumstellar gas, the dust dynamics are more simple. If 
(3 for a given dust size is greater than 0.5, then dust is rapidly removed by radiation 
pressure as in Case 1. If /3 is < 0.5 then Poynting- Robertson drag is the dominant 
removal mechanism. F or /3„ ~ 1 g cm~^ and < Qnhs > 1, the P-R drag timescale can 



be parameterized (e.g. iBackman and Parescdll993l : lAugereau et al.lll999l ): 



where s is the grain size. 

To derive the drag timescale for grains, we consider typical parameters for tau (the 
grain's distance from the star), < Qabs > (the absorption coefficient), and ps (grain 
volume density). Most sources have weak or negligible 8 fim excesses an d stronger 24 



fj,m ex cesses; thus, their dust temperatures are probably < 125-250 K (ICurrie et al. 



2008al ). While 120 K dust can produce a strong 24 fim excess but weak/no 8 fim 
excess, colder dust with 24 fim emission further down the Wien tail can also produce 
these trends. Given the lack of constraints on dust's temperature and location, we 
calculate tp^R for a range of distances - 1, 10, and 30 AU - which cover the range 
of distances for dust at 100-200 K for most stars listed in Table |6] assuming simple 
radiative equilibrium. 

As shown in Table [6], P-R drag removes small grains surrounding early-type IC 348 on 
much timescales shorter than the ~ 2-3 Myr age of the cluster. For 1 fim grains, the 
removal timescales range from being several hundred years to ~ 100,000 years for IDs 
6,8, 20, and 30. In all cases, the removal timescale is less than ^ l/20th of the median 
age for IC 348 stars. Thus, without a replenishment source, dust grains surrounding 
these stars are removed rapidly, on timescales much less than the age of the star. 

The effect of P-R drag on dust surrounding late-type stars is less clear. If the dust or- 
bits closely to the stars (1 AU), then P-R drag can remove the dust on rapid timescales. 
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However, for both ID-39 and ID-1939, it is not clear that the dust requires a replen- 
ishment source if it is located at much larger distances (30 AU). For ID-1939, dust at 
intermediate di stances (10 AU) could also be remnant protoplanetary dust which are 



plausible. The iRobitaille et all ((20061) models that best fit the IDs 39 and 1939 SED 
assumes inner hole sizes of ~ 100 AU and 31 AU, respectively. Thus, the dust need 
not be second generation. 

In summary, the dust grains surrounding the AFG stars listed in Table must be 2nd 
generation. This result is independent of whether the disk has no gas or has residual gas. 
Because the dust must be 2nd generation, disks surrounding these stars must be debris disks. 
Other characteristics of these disks - lack of gas accretion and minimum required dust mass - 
support this conclusion. The only situation in which the dust can be primordial is if the disk 
retains an enormous reservoir of gas, >> 100-1000 M^. In such a case, the gas-dust coupling 
is strong enough that the removal timescale from radiation pressure may be comparable 1-3 
Myr (Takeuchi, pvt. comm.). However, the extreme gas to dust ratios implied by these 
conditions, >> 10^-10^, would be highly unusual if not unphysical. Spectroscopy of sources 
listed in Table E] fail to identify any with active gas accretion, which would be expected if 
the disks had an enormous reservoir of gas left. By far the most simple explanation is that 
these sources are very gas poor or gas free. The dust responsible for their MIPS excesses 
then most likely requires a replenishment source: collisions between larger bodies. Thus, 
these disks fit the standard definition of a debris disk. 

In contrast, it is not clear if dust surrounding late type stars (e.g. ID-39 and ID-1939) 
must be second generation. Both ID-39 and ID-1939 lack evidence for accretion and have 
a low disk mass. However, for ID-39 it is not clear that its dust must second generation if 
the dust is located far from the star and is large {(3 < 0.5). The dust surrounding ID-1939 
need not be second generation if it is located beyond ~ 1 AU of the star and or is entrained 
in a substantial reservoir of gas (> 1 Mq). Both of these situations ar e at least plausible in 



princ iple. Other removal mechanisms such as corpuscular wind drag (jPlavchan et al.l 12005 



20091 ) may be important for late-type, low-mass stars are neglected here. Analyzing dust 
removal timescales after including these effects, may affect our conclusion that the dust 
surrounding late-type stars need not be from a debris disk. 

Based on the above analysis, we classify all AFG stars with characteristics typical of 
those listed in Table [6] (mid IR colors, lack of accretion signatures) as debris disk candidates. 
Debris disk candidates surrounding K and M stars with similar characteristics could be either 
debris disks or evolved primordial disks/transitional disks with very low dust masses. To be 
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conservative, we classify them as homologously depleted disk: 



4. Global Analysis of the IC 348 Disk Population 

We now combine our new optical spectroscopy and MIPS photometry of IC 348 stars 
with archival data to investigate the properties of the entire IC 348 disk population. Our 
analysis has several goals. To investigate the frequency of gas accretion as a function of 
spectral type, we combine our analysis of optical spectra with archival spectroscopic data. 
We consider whether the frequency of signatures of circumstellar gas depends on the spectral 
type of stars, a trend consistent with a stellar-mass dependent dispersal of gas disks. 

Second, we investigate how the rate of dust evolution varies with stellar mass. For a 
simple, empirical probe of disk evolution, we first analyze the Ks-[8] and Ks-[24] colors as 
a function of spectral type and IRAC slope. We consider the evidence for a spectral-type 
dependent dust luminosity at 8 /xm and 24 /im. To further see whether disks surrounding 
stars of different stellar masses evolve at different rates, we conclude by constraining the 
relative frequency of disks in different evolutionary stateqj. 



4.1. Gas Disk Evolution as a Function of Stellar Mass 



To investigate gas disk evolution in IC 3 48, we analyze data present ed in this work, 
archival data summarized in Lada et al. ( 2006 ). and data from Dahm ( 20081). We use empir- 
i cal cr iteria to identify accretors (primarily the Hq, equivalent width) from Iwhite and Basri 
(120031 ) for stars later than GO. Stars with very early spectral types (BAF) have a deep Hq, 
absorption hues produced fro m the stellar photosphere. For these stars, we use our results 
and those from iDahml (120081 ) to identify stars with Hq, emission line reve rsals consistent 
with accretion. Da t a for late- type star s summ arized in iLada et al.l (120061 ) typically draw 
from lLuhman et al.l (120031 ) and iLuhmanI (119991 ). Some sources lacking Hq, equivalent widths 
from optical spectroscopy have ne ar-IR spectra that either confirm their s tatus as accr etors 
(ILuhman et al.l l2005l : iDahml |2008| ) or rule out clear evidence of accretion (IDahml |2008| ) . We 
include these sources in our statistics. We remove from the sample all other stars lacking 



•^While the disk models best fitting the SEDs have large inner holes, the low disk v s. photosphere contra st 
afforded by K and M stars means that the disk mass may simply be extremely low ( Ercolano et al. 2009[ ). 

''We do not include sources that lack ncar-IR photometry because SED modeling of such sources (e.g., 
IDs 746, 2096) was unsuccessful in constraining their evolutionary states 
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published Hq, equivalent widths. 



Figure [TT] (left panel) plots the Hq, equivalent width versus spectral type. Stars identified 
as accretors are shown as black dots while those lacking evidence for accretion are shown 
as grey dots. Stars without Hq, in emission that lack evidence for accretion are shown 
along a line corresponding to EW(Hq)= -15 A(in absorption). For reference, we draw lines 
corresponding to EW(Hq)=0 (the division between emission and absorption) and 10 A, which 
is the division typically used to distinguish between classical T Tauri stars and weak-line T 
Tauri stars. According to our criteria, several stars with spectral types earlier than K5 have 
EW(Hq) < 10 a but are identified as accretors, while others with spectral types later than 
M3 have EW(Hq) > 10 A but are not accreting. 

Figure [H] (right panel) plots the frequency of accretors as a function of spectral type. 
We divide the sample into three spectral type bins: sources earlier than K3, between K3 
and MO, and between MO and M2.5. For a 2-3 Myr-old clust er, these divisions c orrespond 
to stars with masses > 1.4 Mq, 0.8-1.4 Mq, and 0.5-0.8 M© JSaraffe et al.lll998h . Sources 
with probable masses > 1.4 M© have the lowest frequency of accretion (5/29: 17.2+^_4.8%). 
Nearly half of the stars with probable masses of 0.8-1.4 Mq and 0.5-0.8 Mq show evidence 
for gas accretion (7/17, 41.2+i2.2_^^^%^ 43.9+^-«_7.3%, respectively). 

In deriving these frequencies, we acknowledge that usin g Hq as a diagnostic of accretion, 
even following the prescription of IWhite and Basril (120031 ). is prone to some uncertainties. 
The contrast between the Hq emission line and the photosphere is stronger for later type, 
lower- luminosity stars. The large number of stars later than ~ M5 lacking optical spec- 
tra yields an uncertain frequency of accretors for the lowest-mass stars. Nevertheless, the 
observed trend of accretion fr equency is consist ent with the trend of infrared excess as a 
function of spectral type from iLada et al.l (l2006l ) at a ~ 1-2 a significance: the frequency 
is lowest for early type, high-mass stars and highest for 0.6-1 Mq stars. The results are 
consistent with a stellar-na ass dependent dispersal of gaseous circumstellar disks (see also 
Kennedy and KenyonI |2009| ) . 



4.2. Circumstellar Dust Evolution as a Function of Stellar Mass 

To investigate the dust luminosity as a function of stellar mass, we compare the intrinsic 
8 fim and 24 /zm luminosity for each star relative to its intrinsic K-band luminosity. The 
dereddened K5-[8] colors as a function of spectral type and IRAC slope are shown in the 
left panel of Figure [121 The figure includes photometric errors for each source along with 
the locus for photospheric colors from the Kurucz-Lejeuene stellar atmosphere models. The 
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diskless, anemic, and thick sources are clearly separated into a blue to red sequence of colors 
across the entire range of spectral types. Diskless sources (black dots) have Ks-[8] colors 
ranging from ~ to 0.2 from AO to M6/M7. Anemic sources (grey dots) are slightly redder 
(Ks-[8] ~ 0.2-1.25) and have weak excesses. Finally, thick sources (light grey dots) are 
systematically redder, with Ks-[8] ~ 1-3. 

The dereddened K^-fS] color qualitatively shows t he same sep a ration into thick, ane- 
mic, and d i skless classes as the observed a[3.6]-[8] from iLada et al.l (120061 ). Together with 
Lada et al.l (120061 ) . our results show that the median 8 iim excess is larger for later-type 
stars. This result suggests that inner disks around early type, high-mass stars evolve faster 
than around solar and subsolar-mass stars. 

The distribution of dereddened Ks-[24] colors (right panel) is more complicated. There 
is considerable overlap in 24 iim excesses between sources with varying levels of 8 iim excess, 
especially for stars with spectral types between Ml and M7. The 'thick' and 'anemic' stars 
in this spectral type range have a wide range (2-6 mag) of 24 /xm excesses. 

The distribution of Ks-[24] colors suggests that disks around later- type stars are more 
luminous relative to the star even tho ugh the contrast betw een the star and a disk of a given 
dust mass is likely lower for M stars (lErcolano et al.ll2009l ). The fractions of MlPS-detected 
stars with strong 24 /xm excesses (Ks-[24] > 4) for stars earlier than K3 (> 1.4 M©), between 
K3 and MO (0.8-1.4 Mq), and between MO and M2.5 (0.5-0.8 Mq) are 5/18 (27.8%), 7/10 
(70%), and 20/24 (83.3%). 

Despite of the large number of stars lacking detections, this general trend clearly holds 
for stars with masses > 0.5 M©. The 24 /im upper limits for the 12 stars earlier than K3 
without MIPS detections constrain the stars' excesses to be < 3-4 magnitudes. Assuming 
that all stars between K3 and MO currently lacking MIPS detections (10) have weak excesses 
or no excesses, more than 35% of these stars must have strong excesses (7/20), a percentage 
that is still much higher than that for earlier type, higher-mass stars (4/31; 12.9%). 

To compare the median level of excesses for high/intermediate to low- mass IC 348 stars, 
including sources with 24 iim upp e r limits, we use the ASURV Survival Analysis package 
Rev 1.2 (llsobe and Feigelsonl Il986l : iLaValley et al.lll992l ). Because our goal is to compare 
the Kg- [24] colors for two stellar mass bins, we treat the data as univariate and compare 
the two populations using the Gehan's General Wilcoxon Test (Permutated Variance and 
Hypergeometric Variance), the Peto-Prentice Generalized Wilcoxon Test, and the Kaplan- 
Meier Estimator. The three Wilcoxon tests yield the probability that the high/intermediate 
and low-mass populations are drawn from the same population. The Kaplan- Meier Estimator 
yields the likely median value of the Ks-[24] colors for both populations while accounting 
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for non detections. We divide the sample into stars with spectral types earlier than K5 and 
those with spectral types between K5 and M2.5. For a 2.5 Myr-old cl uster, these spectra l 



type ranges correspond to stars with > 1 Mq and = 0.5-1 Mq (IBaraffe et al.lll998l ). 



The ASURV statistical tests indicate that a) the excesses for high/intermedate mass and 
low-mass stars have statistically significant differences and b) the excesses are stronger for 
low-mass stars. The Wilcoxon tests find that the two populations have a probability of less 
than 0.14% (0.001-0.0014) of being drawn from the same parent population. The Kaplan- 
Meier estimator yields a small median Ks-[24] color for high/intermediate- mass sources (K^ 
-[24] ~ 1.55 ± 0.31); the median color for low-mass sources is ~ 3.19 ± 0.40. Very similar 
results are obtained if the division between samples is made anywhere from GO to K5. Thus, 
the late type, low-mass IC 348 population has larger 24 fim excesses from disks. Under 
the general picture that the luminosity of disks at a given wavelength decline with time, 
these results imply that disks around low-mass stars are less evolved than those around 
high/intermediate- mass stars. 

For this analysis to be completely reliable, the censoring of the data must be 'random'. If 
nondetections are primarily due to high nebular background emission at the source positions 
then the censoring is random. If the sources lacking detections are intrinsically fainter, the 
censoring is not random About 56% and 51% of the stars with spectral types earlier than K5 
and between K5 and M2.5 are detected at 24 /im. Because these frequencies are comparable, 
the censoring of data is very close to being random for stars earlier than M2.5. The frequency 
of 24 fim detections drops substantially for stars later than M2.5 (~ 30%). Comparing the 
disk luminosities for > 0.5 M© stars with those for even lower- mass stars is then impeded 
by sample incompleteness for the latter population. 



4.3. Evolutionary State of Disks vs. Stellar Mass in IC 348 

We now consider the frequencies of disks in different evolutionary states as a function 
of stellar mass. For all stars, we quantify the evidence for disk emission at 24 /im. Figure 
|2] suggests that stars with K^- [24] colors > la{Ks, [24]) away from the distribution of pho- 
tospheric sources have significant excesses. Sources with smaller dereddened K<j-[24] colors 
are then identified as diskless photospheres. To guide our identification of primordial disks, 
transitional disks, and debris disks, we use the results from §3.3.3 and §3.4. For stars later 
than KO, we use the semi-empirical division based on Figure [9] to identify primordial disks 
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and transitional disk^. For stars earlier than KO, we use the results of §3.3.3 and §3.4 to 
select primordial disks, transitional disks, and debris disks. 

Figure [T3] shows the frequency of primordial disks, transitional disks, and debris disks 
as a function of spectral type/stellar mass. We divide the sample into the same three stellar 
mass bins used to compute the gas accretion frequency: > 1.4 M© (earlier than K3), 0.8-1.4 
M© (K3-M0), and 0.5-0.8 Mq (M0-M2.5). We do not compute the relative disk fractions 
for later type, lower-mass stars where completeness at 24 fim drops to ~ 30%. 

Based on our analysis, the disk evolutionary states in IC 348 depend on stellar spectral 
type. Half of the MlPS-detected disks around intermediate/high-mass stars are likely debris 
disks (dashed line/diamonds). Primordial disks (solid line/crosses) comprise about 15% 
of the total disk population. The MlPS-detected disk populations for stars with probable 
masses between 0.5 Mq and 1.4 Mq show very different distributions. Primordial disks 
comprise greater than 70% of the MlPS-detected disk population for stars in this mass 
range. The fact that stars in this mass range have (relative to the stellar photosphere) 
more luminous disks than their high/intermediate- mass counterparts is consistent with our 
determination that these stars typically have primordial disks. 

The frequency of transitional disks (dotted line) for the MlPS-detected disk popula- 
tion ranges from ~ 15% for 0.5-0.8 M© stars, to 20% for 0.8-1.4 M© stars, to 35% for 
high/intermediate- mass stars. This frequency is slightly higher than that computed for 
Taurus- Auriga (< 10%). Moreover, the transition disk frequency exceeds the frequency 
of primordial disks for high/intermediate-mass st ars, a result simila r to that found for 



solar /subsolar- mass stars in 5 Myr-old NGC 2362 (jCurrie et al.ll2009l ). Depending on the 
nature of IC 348 sources currently lacking MIPS detections, mid-IR colors presented in Fig- 
ure [12] indicate that the transition disk population for the lowest-mass stars (M^ < 0.5 Mq) 
may also comprise a significant fraction of the total disk population (> 30%). 



5. Summary and Discussion 

5.1. Summary of Results 

Using new optical spectra and new, deep MIPS 24 fim and 70 fim photometry we 
investigated the disk population of the 2.5 Myr-old IC 348 Nebula. Combining these data 



^Within the transitional disk category, we include sources first identified in §3.3.3 as debris disk candidates 
from Figure [9] because the dust around these sources need not be 2nd generation (see §3.4.3). Separating 
these debris disk candidates from transitional disks does not qualitatively change our results. 
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with previous work from iLada et al.l (120061 ) . we analyzed optical spectroscopic data for all 
stars, performed SED modeling, and examined the mid-IR colors of MlPS-detected members 
to constrain the disks' evolutionary states and probe how disk properties vary with stellar 
mass. Our study yields the following major results: 



IC 348 sources with MIPS-70 fim detections include fiat-spectrum protostars and 
pre-main sequence stars with optically-thick, luminous far IR disk emission. Some 
optically-thick disks show evidence for depleted inner disks; others do not. 

IC 348 stars with new and previous MIPS detections have disks in evolutio nary states 
rangi r ig from primordial disks t o debris disks. In a greement with recent work (ILada et al, 



20061 : [Hernandez et al. 1 120071 : ICurrie et al.l l2009l ) , we find evidence for two separate 



pathways from primordial disks to debris disks. In homologously depleted disks, disks 
deplete their reservoir of small dust grains at all disk locations simultaneously. In a 
second sequence ('disks with inner holes'), disks clear their supply of small dust grains 
from the inside out. 

At a ~ 1-2 a significance, signatures of circumstellar gas accretion are more frequent 
for solar and subsolar-mass stars than for high/intermediate-mass stars. This result is 
consistent with a gas disk dispersal timescale that is shortest for high/intermediate- 
mass stars. 

The mid-IR disk luminosities of MlPS-detected disks are stellar-mass dependent. Rela- 
tive to the stellar photosphere, 24 yum emission from disks is lower for high/intermediate- 
mass stars than for solar /subsolar-mass stars. If disks generally decline in luminosity 
as a function of time, this result implies that disks around high/intermediate- mass 
stars evolve faster. 

The evolutionary states of MlPS-detected disks are also stellar mass dependent. Most 
MlPS-detected disks surrounding high/intermediate- mass stars stars appear to be de- 
bris disks; primordial disks comprise only ~ 15% of the disk population for M^ > 1.4 
Mq. In contrast, most MlPS-detected disks around solar and subsolar-mass stars are 
primordial disks. For stars of all masses, transitional disks (homologously depleted or 
disks with inner holes) comprise ~ 15-35% of the MlPS-detected disk population. 



5.2. The Evolutionary State of Anemic Disks 



Our SED modeling and investi gation of the mid -IR colors of IC 348 stars clarifies the 
nature of anemic disks identified by iLada et al.l (120061 ) . Among late type stars, most anemic 



-29- 



disks are homologously depleted disks or are disks with inner holes. Several anemic disks 
are probably primordial disks. Among earlier- type stars with probable masses ^1.4 Mq, 
anemic disks comprise a broader range of evolutionary states. At least two anemic disks 
(surrounding IDs 6 and 8) are probably debris disks. In contrast, ID-31 has strong mid-IR 
emission and gas accretion signatures more consistent with a transitional disk, specifically 
one with an inner hole. 



Thus, as pointed out by ICurrie et al.l (120091 ) . a single IRAC slope serves as a useful 



first-order probe of disk properties, but a full analysis of IRAC and MIPS data combined 
with gas accretion signatures is required for an accurate taxonomy of disks. Our results 
indicate that a full analysis is especially necessary if the stellar population includes both 
high/intermediate-mass stars and lower-mass stars. Given the diversity in disk properties 
for 2.5 Myr-old IC 348, we suggest that a single fiux slope may be best suited as a probe of 
disk evolution in the youngest regions (e.g. Taurus; NGC 1333), where less diversity in disk 
states (e.g. few if any debris disks) is expected. 



5.3. Constraints on the Primordial-to-Debris Disk Transition 



Recent Spitzer studies have placed strong constraints on the timescales for the evolution 
of primordial disks into debris disks. By 5 Myr, most high/intermediate-mass stars and 
solar- mass stars either lack evidence for a disk or h a ve disk properties suggestive of a debris 
disk (jCarpenter et al.l l2006l : iHernandez et al.ll2008l : ICurrie et al.l 120091 ). At this age, many 
(detectable) subs olar-mass stars have weaker levels of rnid-IR disk emission than typical 
primordial disks (ICurrie et al.l l2009l : iDahm and Carpenten |2009|) . 



This work shows that by 2.5 Myr most disks around high/intermediate-mass stars are 
either actively leaving the primordial disk phase or have already reached the debris disk 



phase. Together w ith previous results (e.g. ICarpenter et al.ll2006l : IHernandez et al. 112007 



Currie et al.l 120091 ) . our results clearly support a stellar- mass dependent timescale for the 
disappearance of primordial disks and the emergence of debris disks. Within the context 
of planet formation, gas giant planets around high/intermediate- mass stars have much less 



time to form than around low-mass stars. 



high/intermediate- mass stars (jjohnson et al. 



f gas giant planets are more frequent around 



20071). their formation mus t be very rapid and 



efficient. Only recently have realistic models (IKenyon and Bromleyll2009t ) that form gas giant 



planets via core accretion been successful in forming the cores of such planets in the short 
timescales implied from this paper (2.5 Myr). Some trends in exoplanet properties, such as 
the semimajor axis distribution, may emerge from the competin g effects of core formation 
timescales and primordial disk dispersal timescales (j Currie 2009). 
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The fraction of IC 348 disks in a transitional phase serves as a useful contrast to results 
obtained for younger clusters like Taurus and older clusters like NGC 2362. Based primarily 
on I RAS data, the computed frequency of transitional disks in Taurus is small, << 10% 



[e.g. 



Skrutskie et al.l Il990l : ISimon and Pratd Il995l : IWolk and Walter! Il996l ) . These authors 
then argued that the lifetime of transitional disks is ~ 1-10% the age of Taurus: ~ 0.01-0.1 
Myr. If this inference were correct, the total disk lifetime may be several Myr, but most of 
the disk is dissipated rapidly. 

However, analysi s of Spitzer data fo r 5 Myr-old NGC 2362 finds a much higher frequency 
of transitional disks (ICurrie et al.ll2009l ). Both types of transitional disks (inner holes, ho- 
mologously depleted) greatly outnumber primordial disks. This result argues for far longer 
typical transition timescale (~ 1 Myr). 

In the intermediate age IC 348, the frequency of transition disks is intermediate between 
the frequency for 1 Myr-old Taurus and 5 Myr-old NGC 2362 for solar/subsolar-mass stars. 
Thus, the number of transition disks relative to primordial disks around solar/subsolar-mass 
stars appears to be an increasing function of stellar age. This result is expected if the typical 
transition timescale is an appreciable fraction of the typical primordial disk lifetime. Thus, 
the transition disks in Taurus and IC 348 could remain in such a state for an extended 
period. However, it is possible that transition disk lifetimes, like primordial disk lifetimes, 
also have an intrinsic dispersion (i.e., a gaussian distribution of lifetimes), where sources in 
Taurus represent those that make the primordial-to-debris disk transition fastest and sources 
in 5-10 Myr-old clusters make the transition the slowest. Quantifying the range of times 
a disk spends in a transitional phase and whether the typical timescale depends on stellar 
mass requires observations of many more 1-10 Myr-old clusters. 

While debris disks typically only have excess emission longwards of ~ 20 /xm, the debris 
disk population in IC 348 contains at least two sources (ID-6 and ID-8) that have warmer dust 
more indicative of terrestrial planet formation than icy planet formation. Thus, warm debris 
disks consistent with the observable signatures of terrestrial planet formation may emerge as 
early as ~ 2.5 Myr around high/intermediate-mass stars. Observations of more 2-5 Myr-old 
clusters are needed to determine if the frequency of warm debris disks is higher than the ~ 
4% derived for 10-15 Myr-old clusters (ICurrie et al.ll2007al : ICurriell2008l ). The existence of 
both warm debris disks and colder debris disks (lacking IRAC excess emission) suggests that 
even at 2.5 Myr debris disk populations may exhibit a range of dust temperatures consistent 
with a range of locations over which the debris-producing stages of planet formation are 
ongoing. 
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Table 1. Spectroscopic data for Selected IC 348 sources 



ID Number Spectral IRAC SNR EW 


Emission? 


FWHM EW Ka Rc(est.) R* 




Flux(Ha) 


M/MlD-21 


Type 


Slope 


Ha (A) 




H<,(A) 


Em. (A) 


(R© 


) (M© 


) (W m-2) 




2 A2 


thick 143 


-10.7 


maybe 




< 0.015 7.2 


4.5 


3.1 


< 3.8x10-1'' 


<0.07 


3 AO 


diskless 170 


-6.3 


no 




< 0.02 7.3 


4.0 


3.0 


< 4.6x10-1^ 


<0.08 


6 G3 


anemic 105 


-1.0 


probably not 




< 0.1 9.0 


3.9 


2.3 


< 4.8x10-1^ 


< 0.10 


8 A2 


anemic 156 


-10.9 


no 




<0.02 8.51 


2.4 


2.3 


<1.5xl0-" 


<0.02 


19 A2 


thick 132 


-10.0 


yes 


2.9 


0.64 9.15 


1.7 


2.1 


2.6x10-1'^ 


0.3 


21 KO 


anemic 80 


5.0 


yes 


7.8 


5.0 10.48 


2.5 


1.8 


6x10-1'^ 


1 


30 FO 


diskless 123 


-6.2 


no 




<0.02 10.06 


1.6 


1.6 


<3.5xl0-i** 


<0.004 



Note. — Negative EW{Hq) indicate absorption and positive values indicate emission. IC348 ID-19 has a central emission 
core in the deep absorption line. The Ha line for ID-6 is weak for G3 stars, though there is no clear central emission feature 
and the star is very x-ray luminous. 



Table 2. IC 348 Stars with New MIPS 24 fim Detections 



ID RA DEC OffsetC) Spectral Type IRAC Slope [24] cr[24] Ks-[24] (dereddened) 



4 


56, 


.1300 


32, 


,1061 


1.80 


FO 


3 


6.402 


0.097 1.24 


39 


56, 


.2573 


32 


.2411 


3.68 


K4 


3 


8.818 


0.225 0.87 


47 


55, 


.9813 


32, 


,1590 


1.43 


KO 


3 


10.488 


0.280 -0.42 


49 


55, 


.9900 


32, 


,0271 


3.72 


M0.5 


3 


8.076 


0.430 0.53 


53 


56, 


.0684 


32, 


.1653 


0.29 


KO 


3 


9.545 


0.257 0.62 


56 


56, 


,0208 


32, 


,1649 


1.62 


K3.5 


3 


9.838 


0.195 0.45 


83 


56, 


,1559 


32 


.1503 


0.18 


Ml 


1 


5.733 


0.286 5.05 


85 


56, 


,1172 


32, 


.2668 


0.45 


M3.25 


3 


9.078 


0.350 1.84 


88 


56, 


.1365 


32, 


.1544 


2.97 


M3.25 


2 


6.819 


0.425 4.00 


91 


56, 


.1634 


32, 


,1624 


1.90 


M2 


1 


8.264 


0.265 2.63 


95 


56, 


.0913 


32, 


.2032 


1.16 


M4 


3 


9.118 


0.311 1.88 


124 


55, 


.9776 


32, 


.0084 




M4.25 


3 


9.144 


0.190 2.19 


167 


56, 


.1716 


32, 


,1695 


0.47 


M3 


2 


7.246 


0.565 3.93 


168 


56, 


,1306 


32 


.1797 


1.15 


M4.25 


1 


7.145 


0.263 4.37 


193 


56, 


.1584 


32, 


,1936 


0.85 


M4 


1 


7.141 


0.182 4.76 


261 


55, 


.9526 


32, 


,2308 


1.96 


M5 


2 


10.037 


0.203 2.65 


300 


56, 


.1624 


32, 


,0555 


1.49 


M5 


1 


7.735 


0.265 5.23 


336 


56, 


,1349 


32, 


.0576 


0.95 


M5.5 


1 


8.058 


0.335 5.27 


341 


56, 


.0541 


32, 


.2210 


1.08 


M5.25 


1 


8.853 


0.304 4.28 


365 


56, 


,0426 


32, 


,1263 


1.26 


M5.75 


2 


9.383 


0.340 4.01 


373 


56, 


.1166 


32, 


,0888 




M5.5 


2 


8.418 


0.107 5.31 


402 


56, 


.1898 


32, 


.3055 


0.41 


M5.5 


2 


8.611 


0.174 4.92 


407 


56, 


.2672 


32, 


.0846 


0.25 


M7 


1 


9.578 


0.164 4.89 


746 


56, 


.2082 


32, 


.1041 


1.40 


M5 


1 


7.683 


0.260 ? 


1124 56, 


.2364 


32, 


,2844 


0.86 


M5 


2 


8.406 


0.198 3.53 


1939 56, 


.2198 


32, 


,0158 


1.11 


M4.75 


3 


10.584 


0.313 0.82 


2096 56, 


.0539 


32, 


,2234 


0.39 


M6 


1 


8.713 


0.233 ? 



Note. — The IRAC sl ope follows the conv ention of iLada et al.l ll2006h : 1 = 'thick', 2 = 'anemic', and 3 = 'diskless'. The 
equivalent designations in lCurrie et al. 1 1 I2OO9I) are the following: 1 = 'strong IRAC, 2 = 'weak IRAC, and 3 = 'photosphere'. 
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Table 3. IC 348 Stars With MIPS 70 /xm Detections 



IC348-ID Spectral Type Source Type F70 (mjy) (t(F7o) (mjy) OffsetC) 



13 


M0.5 


1 


1408 


103.7 


2.85 


31 


Gl 


2 


556.2 


140.3 


2.86 


51 




1 


3972.0 


127.0 


1.74 


67 


M0.75 


2 


125.5 


18.4 


2.08 


10343 


M3.75 


1 


115.6 


15.9 


1.96 


10352 


Ml 


1 


322.1 


25.1 


2.37 



Table 4. Catalog of IC 348 Stars 



ID RA DEC J H Ks [3.6] cr[3.6] [4.5] o-[4.5] [5.8] cr[5.8] [8] a[8] [24] o-[24] [70] o-[70] ST IRAC Slope 



1 56.1425 32.1630 6.790 6.640 6.510 6.740 0.010 6.540 0.020 

2 56.1474 32.1679 7.950 7.530 7.250 7.090 0.020 6.810 0.020 

3 56.2110 32.3185 8.370 7.900 7.660 7.530 0.010 7.610 0.030 

4 56.1300 32.1061 8.340 8.020 7.860 7.790 0.010 7.740 0.020 

5 56.1084 32.0751 10.070 8.870 8.140 6.970 0.010 6.520 0.010 



6.580 0.020 6.500 0.030 1.874 -9 -2.267 -9 15 3.10 3 

6.460 0.040 5.820 0.040 3.375 0.038 -2.176 -9 22 3.20 1 

7.470 0.040 7.500 0.040 7.086 0.141 -1.009 -9 20 3.90 3 

7.660 0.020 7.730 0.030 6.402 0.097 -0.898 -9 30 2.30 3 

6.320 0.020 5.630 0.020 2.655 0.012 -0.699 -9 48 7.70 1 



Note. — Full catalog of IC 348 sources from flTada et aL I l l2006h including MIPS photometry for newly-detected stars from 
this work. The column labeled "ST" indicates the numerical spectral type (e.g., 15=B5, 22=A2). 



Table 5. Evolutionary State of Disks with New 24 fim Detections 



ID Spectral Type IRAC Slope Ks-[8] (dereddened) Ks-[24] (dereddened) Disk State 



4 


FO 


3 




-0.01 




1.24 






4 


39 


K4 


3 




0.08 




0.87 






4 


47 


KO 


3 




0.10 




-0.42 






5 


49 


MO. 5 


3 




0.24 




0.53 






5 


53 


KO 


3 




0.01 




0.62 






4 


56 


K3.5 


3 




0.02 




0.45 






4 


83 


Ml 


1 




2.00 




5.05 






1 


85 


M3.25 


3 




0.26 




1.84 






2 


88 


M3.25 


2 




0.60 




4.00 






3 


91 


M2 


1 




1.71 




2.63 






2 


95 


M4 


3 




0.39 




1.88 






2 


124 


M4.25 


3 




0.55 




2.19 






2 


167 


M3 


2 




1.16 




3.93 






2 


168 


M4.25 


1 




1.91 




4.37 






1 


193 


M4 


1 




1.76 




4.76 






1 


261 


M5 


2 




0.66 




2.65 






2 


300 


M5 


1 




1.97 




5.23 






1 


336 


M5.5 


1 




1.99 




5.27 






1 


341 


M5.25 


1 




1.40 




4.28 






1 


365 


M5.75 


2 




1.55 




4.01 






1 


373 


M5.5 


2 








5.31 






3 


402 


M5.5 


2 




1.32 




4.92 






1 


407 


M7 


1 




2.15 




4.89 






1 


746 


M5 


1 
















1124 M5 


2 




0.70 




3.53 






2 


1939 M4.75 


3 




0.39 




0.82 






4 


2096 M6 


1 
















Note. — The IRAC slope follows the convention oflLada 


et all 


|2006|): 1 


= 'thick', 2= 'anemic', and 3 = 'diskless'. The disk 


states are 1 = primordial, 2 = 


= homologously depleted, 3 = disks with inner 


holes, 4 = debris disk candidates, 5 = photospheres 


(no 


disk). 




















Table 6. 


Dust Removal Timescales and Dust Masses for Representative Debris Disk 














Candidates 




ID 


IRAC" Spectral log(L4/LQ 


) M/Mq tp^ (Myr) 








Disk'' 






Slope Type 






(1, 10, 30 AU) 






Model 


(xIO^^^Mq) 


6 


Anemic G3 


1.209 


2.6 


4.3x10^^, 


4.3x10" 


0.039 7.1 


3006897-5- 


■ 9 8.77 


8 


Anemic A2 


1.509 


2.8 


2.2x10"^, 


2.2x10" 


^, 0.020 13.1 


3004289-5- 


■ 9 2.48 


20 


Diskless Gl 


0.719 


2.6 


1.3x10"'*, 


0.013, 0.120 


2.3 


3001439-5- 


9 2.09 


30 


Diskless FO 


0.806 


2.6 


l.lxlO-*, 


0.011, 0.098 


2.8 


3003786-5- 


9 2.95 


39 


Diskless K4 


0.229 


1.3 


4.1x10-*, 


0.041, 0.371 


1.5 


3008796-8- 


9 20.0 


1939 Diskless M4.75 


-0.792 


0.2 


4.4x10"^, 


0.44, 3.9 


0.98 


3002293-6- 


■9 14.0 



Note. — Constraints on the dust removal timescales and disk masses for debris disk candidates. The Poynting-Robertson drag timcscalc is 
calculated assuming 1 /j,m-sizcd grains. Disk models listed represent the file number from the iRobitaille et al.| i2006|) grid; dust masses for those 
models are listed as M^.^^^^. a) From fLada et al.| i2006| ) . b) From [Robitaille et al ] 120061) . 
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Fig. 1. — FAST spectra of IC 348 sources. The spectr a have not been smoothed. Each 
source is labled by its catalog number from iLada et aLl (120061 ). its spectral type, and its 
IRAC slope. ID-21 (KO anemic) has a strong, wide Hq emission line; ID-19 (A2 thick) also 
has a Hq, emission core in its absorption line. ID-6 (G3 anemic) has Hq, nearly filled in, 

thniicrVi this is likplv Hiip tn strrmcr phrmTinsr>Vip>rip artivitv Thp ntVipr stars rlparlv hav(=> VI 



-4 -2 2 4 5 8 

Ks-[Z4] 



-2 2 4 6 

Kb-[24] 



Fig. 2.— (left) J-K, vs. K,-[24] for all 451 MIPS-24 fxm detections with 2MASS counterparts. 
The dashed line identifies the locus of photospheric colors for a star reddened by Ay— 2. 
(right) The [24] vs. K5-[24] color-magnitude diagram for the same sources. The long dashed 
line identifies the photospheric locus for a star with Ks-[24]= 0.25: a mid-K star reddened 
by Ay=2 or an unreddened M2-M3 star. The dotted line, short-dashed line, and dash-three 
dots line correspond to the 1 and 2 a limits from the combined and [24] photometric 
uncertainties assuming a ~ 1 magnitude spread in extinction. 
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Fig. 3. — Three color mosaic of the IC 348 Nebula. Blue represents IRAC 8 /im Cycle 1 
data produced by the post-bed pipeline mosaic; green and red colors represent the MIPS-24 
fim and 70 fim emission from the processed data described in this paper. 
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Fig. 4. — Comparisons of our MIPS 24 /im phot ometry/upper limits (left) and photometric 
errors (right) with those from iLada et al.l (120061 ). In the left panel, grey dots identify mag- 
nitudes for sources detected from both datasets and crosses identify the 2a upper limits for 
sources that lack dete ctions in eit h er da taset. The solid line identifies where our [24] pho- 
tometry and that from lLada et al.l (120061 ) would be in perfect agreement in the photon-noise 
limit. The dotted line represe i its th e expected locus of upper limits, given our exposure 
time, based on the iLada et al.l (120061 ) upper limits. In the right panel, grey dots represent 
photometric errors and the solid line identifies where our photometric errors are identical. 




Fig. 5. — [3.6]- [8] vs. [8]- [24] color-color diagram for IC 348 stars with new MIPS 24 [im 
detections. Overplotted are the Icr error bars. For comparison, we show the colors for the 
median Taurus SED (asterisks). (Left) Stars with spectral types earlier than KO are shown 
as dark circles; stars with spectral types later than KO are shown as grey circles. (Right) 
Stars identified as having strong ('thick') IRAC excess emission are dark circles, stars with 
weak ('anemic') IRAC excess emission are grey, and stars lacking IRAC excess emission 
('diskless') are light grey. In bo th plots, a reddening vec tor of Ai,=20 is overplotted, where 
we assume extinctio n laws from llridebetouw et al.l ( 20051) for IRAC and approxim ate the 24 



Flaherty and MuzeroUd (l2007l ) estimate 



[im extinction from iMathid (1l990l ) . We note that 
that the extinction law is much flatter from 8 [im to 24 /im than assumed here. 




Fig. 6. — [3.6]-[24] vs. [24]-[70] color-color diagram for IC 348 stars detected at 70 fim (dots). 
Overplotted are the positions of two rj Cha stars (diamonds), ECH-2 (lower left) and ECH-15 
(upper right). 
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Fig. 7. — SED fits to IC 348 sources detected at 70 fiui using the iRobitaille et al.l (120061 ) 
models. The dotted hue corresponds to the unreddened SE P of the best-fit stella r photo- 
sphere. The solid black line correspond s to the best-fit SEP IRobitaille et al. ( 2006 ) model. 
The light grey lines correspond to other IRobitaille et al. models with a x^-xlest < 10- 
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Fig. 8. — Representative SEDs of IC 348 stars with new MIPS 24 detections. Data 
points are black dots with photometric errors overplotted. The sohd hne corresponds to a 
stellar photosphere approp riate for the star's spectral t ype, the dashed line corresponds to 
the debris disk model from lKenyon and Bromleyl (120041 ). the dotted line corresponds to the 
SED of the median Taurus SED, and the dot-dashed lines correspond to the upper and lower 
quartiles for the median Taurus SED. Each source is identified by its IRAC slope (Strong 
IRAC= 'thick', Weak IRAC= 'anemic', Photosphere= 'diskless), its optical extinction, its 
J2000 position in decimal degrees, and its catalog ID number. 
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Fig. 9. — (top) K5-[8] vs. Ks-[24] color-color diagram for sources having spectral types KO and 
later with new MIPS detections (open symbols) and previously-d etected sources (filled dots) 



later than KO. The colors are dereddened using Ay values from iMuench et al.l (120071 ). For 
sources with new detections, squares correspond to bare photospheres, inverted triangles are 
debris disk candidates, triangles are disks with inner holes, stars are homologously depleted 
disks, and open circles are primordial disks. The disk population is well separated in color- 
color space (dashed and dash-dotted lines). Previously detected sources labeled as 'thick' 
tend to reside in the primordial disk region, sources identified as anemic are located in 
the homologously depleted and inner-holed regions, and 'diskless' sources are located in the 
debris disk/photospheric region, (bottom) The same diagram for sources with spectral types 
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Fig. 10. — SEP of represen tative debris disk candidates. Overplotted are the best-fit SEDs 

(|200(:, soli d line) along with the stellar photosphere(dotted line). Also 



from Robitaille et al 



listed are the 



Lada et al. 



(120061 ) IRAC slope and ID number, the numerical spectral type 
(below the IRAC slope), the J200 source position (in degrees), and the model number for 
the best-fit iRobitaille et al.l (120061 ) model SED. Photometric uncertainties (1 a) are typically 
smaller than the symbol sizes. 



0.0 



Fig. 11. — (left) Hq equivalent width versus spectral type for IC 348 stars. Accretors are 
shown as black dots; nonaccretors are grey dots. Sources without measured EW(Ha) values 
have EW(Hq,) set to -15 A. For reference, we draw a line corresponding to the division between 
Classical T Tauri stars and Weak-hne T Tauri stars, (right) Frequency of accretors as a 
function of spectral type for three spectral type bins. The horizontal error bars correspond 
to the bin size in spectral type - B5-K3, K3-M0, and M0-M2.5 - and the vertical error bar 
identifies the 1 a uncertainty in the frequency of accretors for each spectral type bin based 
on binomial statistics. 
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Fig. 12. — Ks-[8] (left panel) and Ks-[24] (right panel) colors vs. spectral type for sources 
with MIPS 24 fj,m detections. Sources are divided by their IRAC slope: 'thick' sources are 
light grey dots, 'anemic' sources are dark grey dots, and 'diskless' sources are black dots. The 
vertical dotted line divides the sample into stars with probable masses > 1.4 Mq, 0.8-1.4 Mq, 
and < 0.8 Mq. The horizontal solid line shows the locus of photospheric colors according 
to the Kurucz-Lejeune stellar atmosphere models. The long-dashed horizontal line identifies 
sources with a > 20% IR excess. 
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Fig. 13. — Disk evolutionary state from mid-IR colors as a function of stellar mass. The 
frequencies are given for three stellar mass bins (from left to right): MlPS-detected stars 
with probable masses > 1.4 Mq, 0.8-1.4 M©, and 0.5-0.8 M©. 
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Fig. 14. — Atlas of spectral energy distributions of IC 348 stars with new MIPS 24 //m 
detections. Overplotted are the median Taurus SED (dotted line) with upper and lower 
quartiles (dot-dashed line). A terrestrial zone debris disk model is shown as a dashed line with 
mid-IR flux slightly greater than the photosphere (solid line, connected by open squares). 
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